Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes.

BACKGROUND: Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. RESULTS: Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. CONCLUSION: We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Blood DNA methylation profiling identifies cathepsin Z dysregulation in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is characterised by pulmonary vascular remodelling causing premature death from right heart failure. Established DNA variants influence PAH risk, but susceptibility from epigenetic changes is unknown. We addressed this through epigenome-wide association study (EWAS), testing 865,848 CpG sites for association with PAH in 429 individuals with PAH and 1226 controls. Three loci, at Cathepsin Z (CTSZ, cg04917472), Conserved oligomeric Golgi complex 6 (COG6, cg27396197), and Zinc Finger Protein 678 (ZNF678, cg03144189), reached epigenome-wide significance (p < 10-7) and are hypermethylated in PAH, including in individuals with PAH at 1-year follow-up. Of 16 established PAH genes, only cg10976975 in BMP10 shows hypermethylation in PAH. Hypermethylation at CTSZ is associated with decreased blood cathepsin Z mRNA levels. Knockdown of CTSZ expression in human pulmonary artery endothelial cells increases caspase-3/7 activity (p < 10-4). DNA methylation profiles are altered in PAH, exemplified by the pulmonary endothelial function modifier CTSZ, encoding protease cathepsin Z.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Bidirectional Mendelian Randomization and Multiphenotype GWAS Show Causality and Shared Pathophysiology Between Depression and Type 2 Diabetes.

OBJECTIVE: Depression is a common comorbidity of type 2 diabetes. We assessed the causal relationships and shared genetics between them. RESEARCH DESIGN AND METHODS: We applied two-sample, bidirectional Mendelian randomization (MR) to assess causality between type 2 diabetes and depression. We investigated potential mediation using two-step MR. To identify shared genetics, we performed 1) genome-wide association studies (GWAS) separately and 2) multiphenotype GWAS (MP-GWAS) of type 2 diabetes (19,344 case subjects, 463,641 control subjects) and depression using major depressive disorder (MDD) (5,262 case subjects, 86,275 control subjects) and self-reported depressive symptoms (n = 153,079) in the UK Biobank. We analyzed expression quantitative trait locus (eQTL) data from public databases to identify target genes in relevant tissues. RESULTS: MR demonstrated a significant causal effect of depression on type 2 diabetes (odds ratio 1.26 [95% CI 1.11-1.44], P = 5.46 × 10-4) but not in the reverse direction. Mediation analysis indicated that 36.5% (12.4-57.6%, P = 0.0499) of the effect from depression on type 2 diabetes was mediated by BMI. GWAS of type 2 diabetes and depressive symptoms did not identify shared loci. MP-GWAS identified seven shared loci mapped to TCF7L2, CDKAL1, IGF2BP2, SPRY2, CCND2-AS1, IRS1, CDKN2B-AS1. MDD has not brought any significant association in either GWAS or MP-GWAS. Most MP-GWAS loci had an eQTL, including single nucleotide polymorphisms implicating the cell cycle gene CCND2 in pancreatic islets and brain and the insulin signaling gene IRS1 in adipose tissue, suggesting a multitissue and pleiotropic underlying mechanism. CONCLUSIONS: Our results highlight the importance to prevent type 2 diabetes at the onset of depressive symptoms and the need to maintain a healthy weight in the context of its effect on depression and type 2 diabetes comorbidity.

Abdominal obesity is a more important causal risk factor for pancreatic cancer than overall obesity.

Obesity and type 2 diabetes (T2D) are associated with increased risk of pancreatic cancer. Here we assessed the relationship between pancreatic cancer and two distinct measures of obesity, namely total adiposity, using BMI, versus abdominal adiposity, using BMI adjusted waist-to-hip ratio (WHRadjBMI) by utilising polygenic scores (PGS) and Mendelian randomisation (MR) analyses. We constructed z-score weighted PGS for BMI and WHRadjBMI using publicly available data and tested for their association with pancreatic cancer defined in UK biobank (UKBB). Using publicly available summary statistics, we then performed bi-directional MR analyses between the two obesity traits and pancreatic cancer. PGSBMI was significantly (multiple testing-corrected) associated with pancreatic cancer (OR[95%CI] = 1.0804[1.025-1.14], P = 0.0037). The significance of association declined after T2D adjustment (OR[95%CI] = 1.073[1.018-1.13], P = 0.00904). PGSWHRadjBMI association with pancreatic cancer was at the margin of statistical significance (OR[95%CI] = 1.047[0.99-1.104], P = 0.086). T2D adjustment effectively lost any suggestive association of PGSWHRadjBMI with pancreatic cancer (OR[95%CI] = 1.039[0.99-1.097], P = 0.14). MR analyses showed a nominally significant causal effect of WHRadjBMI on pancreatic cancer (OR[95%CI] = 1.00095[1.00011-1.0018], P = 0.027) but not for BMI on pancreatic cancer. Overall, we show that abdominal adiposity measured using WHRadjBMI, may be a more important causal risk factor for pancreatic cancer compared to total adiposity, with T2D being a potential driver of this relationship.

Maternal and perinatal obesity induce bronchial obstruction and pulmonary hypertension via IL-6-FoxO1-axis in later life.

Obesity is a pre-disposing condition for chronic obstructive pulmonary disease, asthma, and pulmonary arterial hypertension. Accumulating evidence suggests that metabolic influences during development can determine chronic lung diseases (CLD). We demonstrate that maternal obesity causes early metabolic disorder in the offspring. Here, interleukin-6 induced bronchial and microvascular smooth muscle cell (SMC) hyperproliferation and increased airway and pulmonary vascular resistance. The key anti-proliferative transcription factor FoxO1 was inactivated via nuclear exclusion. These findings were confirmed using primary SMC treated with interleukin-6 and pharmacological FoxO1 inhibition as well as genetic FoxO1 ablation and constitutive activation. In vivo, we reproduced the structural and functional alterations in offspring of obese dams via the SMC-specific ablation of FoxO1. The reconstitution of FoxO1 using IL-6-deficient mice and pharmacological treatment did not protect against metabolic disorder but prevented SMC hyperproliferation. In human observational studies, childhood obesity was associated with reduced forced expiratory volume in 1 s/forced vital capacity ratio Z-score (used as proxy for lung function) and asthma. We conclude that the interleukin-6-FoxO1 pathway in SMC is a molecular mechanism by which perinatal obesity programs the bronchial and vascular structure and function, thereby driving CLD development. Thus, FoxO1 reconstitution provides a potential therapeutic option for preventing this metabolic programming of CLD.

LongITools: Dynamic longitudinal exposome trajectories in cardiovascular and metabolic noncommunicable diseases.

The current epidemics of cardiovascular and metabolic noncommunicable diseases have emerged alongside dramatic modifications in lifestyle and living environments. These correspond to changes in our "modern" postwar societies globally characterized by rural-to-urban migration, modernization of agricultural practices, and transportation, climate change, and aging. Evidence suggests that these changes are related to each other, although the social and biological mechanisms as well as their interactions have yet to be uncovered. LongITools, as one of the 9 projects included in the European Human Exposome Network, will tackle this environmental health equation linking multidimensional environmental exposures to the occurrence of cardiovascular and metabolic noncommunicable diseases.

Relationship between glucose homeostasis and obesity in early life-a study of Italian children and adolescents.

Epidemic obesity is the most important risk factor for prediabetes and type 2 diabetes (T2D) in youth as it is in adults. Obesity shares pathophysiological mechanisms with T2D and is likely to share part of the genetic background. We aimed to test if weighted genetic risk scores (GRSs) for T2D, fasting glucose (FG) and fasting insulin (FI) predict glycaemic traits and if there is a causal relationship between obesity and impaired glucose metabolism in children and adolescents. Genotyping of 42 SNPs established by genome-wide association studies for T2D, FG and FI was performed in 1660 Italian youths aged between 2 and 19 years. We defined GRS for T2D, FG and FI and tested their effects on glycaemic traits, including FG, FI, indices of insulin resistance/beta cell function and body mass index (BMI). We evaluated causal relationships between obesity and FG/FI using one-sample Mendelian randomization analyses in both directions. GRS-FG was associated with FG (beta = 0.075 mmol/l, SE = 0.011, P = 1.58 × 10-11) and beta cell function (beta = -0.041, SE = 0.0090 P = 5.13 × 10-6). GRS-T2D also demonstrated an association with beta cell function (beta = -0.020, SE = 0.021 P = 0.030). We detected a causal effect of increased BMI on levels of FI in Italian youths (beta = 0.31 ln (pmol/l), 95%CI [0.078, 0.54], P = 0.0085), while there was no effect of FG/FI levels on BMI. Our results demonstrate that the glycaemic and T2D risk genetic variants contribute to higher FG and FI levels and decreased beta cell function in children and adolescents. The causal effects of adiposity on increased insulin resistance are detectable from childhood age.

Supplementation with Iron in Pulmonary Arterial Hypertension. Two Randomized Crossover Trials.

Rationale: Iron deficiency, in the absence of anemia, is common in patients with idiopathic and heritable pulmonary arterial hypertension (PAH) and is associated with a worse clinical outcome. Oral iron absorption may be impeded by elevated circulating hepcidin concentrations. The safety and benefit of parenteral iron replacement in this patient population is unclear. Objectives: To evaluate the safety and efficacy of parenteral iron replacement in PAH. Methods: In two randomized, double-blind, placebo-controlled 12-week crossover studies, 39 patients in Europe received a single infusion of ferric carboxymaltose (Ferinject) (1,000 mg or 15 mg/kg if weight <66.7 kg) or saline as placebo, and 17 patients in China received iron dextran (Cosmofer) (20 mg iron/kg body weight) or saline placebo. All patients had idiopathic or heritable PAH and iron deficiency at entry as defined by a serum ferritin <37 µg/L or iron <10.3 µmol/L or transferrin saturations <16.4%. Results: Both iron treatments were well tolerated and improved iron status. Analyzed separately and combined, there was no effect on any measure of exercise capacity (using cardiopulmonary exercise testing or 6-minute walk test) or cardiopulmonary hemodynamics, as assessed by right heart catheterization, cardiac magnetic resonance, or plasma NT-proBNP (N-terminal-pro hormone brain natriuretic peptide) at 12 weeks. Conclusions: Iron repletion by administration of a slow-release iron preparation as a single infusion to patients with PAH with iron deficiency without overt anemia was well tolerated but provided no significant clinical benefit at 12 weeks. Clinical trial registered with ClinicalTrials.gov (NCT01447628).

Sex-dimorphic genetic effects and novel loci for fasting glucose and insulin variability.

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Expression Quantitative Trait Locus Mapping in Pulmonary Arterial Hypertension.

Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.

Heterogeneous Effects of Calorie Content and Nutritional Components Underlie Dietary Influence on Pancreatic Cancer Susceptibility.

Pancreatic cancer is a rare but fatal form of cancer, the fourth highest in absolute mortality. Known risk factors include obesity, diet, and type 2 diabetes; however, the low incidence rate and interconnection of these factors confound the isolation of individual effects. Here, we use epidemiological analysis of prospective human cohorts and parallel tracking of pancreatic cancer in mice to dissect the effects of obesity, diet, and diabetes on pancreatic cancer. Through longitudinal monitoring and multi-omics analysis in mice, we found distinct effects of protein, sugar, and fat dietary components, with dietary sugars increasing Mad2l1 expression and tumor proliferation. Using epidemiological approaches in humans, we find that dietary sugars give a MAD2L1 genotype-dependent increased susceptibility to pancreatic cancer. The translation of these results to a clinical setting could aid in the identification of the at-risk population for screening and potentially harness dietary modification as a therapeutic measure.

Mendelian randomisation analysis of red cell distribution width in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is a rare disease that leads to premature death from right heart failure. It is strongly associated with elevated red cell distribution width (RDW), a correlate of several iron status biomarkers. High RDW values can signal early-stage iron deficiency or iron deficiency anaemia. This study investigated whether elevated RDW is causally associated with PAH.A two-sample Mendelian randomisation (MR) approach was applied to investigate whether genetic predisposition to higher levels of RDW increases the odds of developing PAH. Primary and secondary MR analyses were performed using all available genome-wide significant RDW variants (n=179) and five genome-wide significant RDW variants that act via systemic iron status, respectively.We confirmed the observed association between RDW and PAH (OR 1.90, 95% CI 1.80-2.01) in a multicentre case-control study (cases n=642, disease controls n=15 889). The primary MR analysis was adequately powered to detect a causal effect (odds ratio) between 1.25 and 1.52 or greater based on estimates reported in the RDW genome-wide association study or from our own data. There was no evidence for a causal association between RDW and PAH in either the primary (ORcausal 1.07, 95% CI 0.92-1.24) or the secondary (ORcausal 1.09, 95% CI 0.77-1.54) MR analysis.The results suggest that at least some of the observed association of RDW with PAH is secondary to disease progression. Results of iron therapeutic trials in PAH should be interpreted with caution, as any improvements observed may not be mechanistically linked to the development of PAH.

Genetic determinants of risk in pulmonary arterial hypertension: international genome-wide association studies and meta-analysis.

BACKGROUND: Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. METHODS: We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses. FINDINGS: A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55-2·08], p=5·13 × 10-15) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42-1·71], p=7·65 × 10-20) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25-1·48], p=1·69 × 10-12; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02-8·05]), despite similar baseline disease severity. INTERPRETATION: This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials. FUNDING: UK NIHR, BHF, UK MRC, Dinosaur Trust, NIH/NHLBI, ERS, EMBO, Wellcome Trust, EU, AHA, ACClinPharm, Netherlands CVRI, Dutch Heart Foundation, Dutch Federation of UMC, Netherlands OHRD and RNAS, German DFG, German BMBF, APH Paris, INSERM, Université Paris-Sud, and French ANR.

Recent advances in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is a rare disorder with a high mortality rate. Treatment options have improved in the last 20 years, but patients still die prematurely of right heart failure. Though rare, it is heterogeneous at the genetic and molecular level, and understanding and exploiting this is key to the development of more effective treatments. BMPR2, encoding bone morphogenetic receptor type 2, is the most commonly affected gene in both familial and non-familial PAH, but rare mutations have been identified in other genes. Transcriptomic, proteomic, and metabolomic studies looking for endophenotypes are under way. There is no shortage of candidate new drug targets for PAH, but the selection and prioritisation of these are challenges for the research community.

Endovascular therapy for vasospasm after aneurysmatic subarachnoid hemorrhage.

INTRODUCTION: Balloon angioplasty and/or selective intra-arterial vasodilator therapies are treatment options in patients with vasospasm after subarachnoid hemorrhage (SAH). We analyzed the effect of balloon angioplasty and/or selective intra-arterial vasodilator therapy in our patients. METHODS: Twenty-six patients (vasodilation group, VDT) were treated with intra-arterial nimodipine. The balloon angioplasty with nimodiopine-group (BAP-N group) comprised 21 patients. The primary endpoint of this study was successful angiographic vessel dilation in vasospastic vessels after balloon angioplasty, together with nimodipine (BAP-N group), compared to intra-arterial vasodilator therapy (VDT group) with nimodipine alone. RESULTS: A significant effect of angioplasty plus nimodipine was found in the central arteries (composite endpoint) with an OR of 2.4 (95% CI: 1.4-4.2], p = 0.002), indicating a chance of improvement of the BAP-N group of more than twice compared to nimodipine infusions alone. Significant advantages for BAP-N-therapy were also encountered in the internal carotid artery (OR 5.4, p < 0.001) and basilar artery (OR 29.7, p = 0.003). A joint analysis of all arteries combined failed to show significant benefit of BAP-N therapy (OR 1.5, p = 0.079), which was also true for cerebral peripheral arteries (OR 0.77, p = 0.367). There was no difference in clinical outcome between both groups. CONCLUSIONS: In SAH patients with vasospasm, a combination therapy of balloon angioplasty and intra-arterial nimodipine resulted in a more than doubled vasodilative effect in the central cerebral arteries compared to the sole infusion of nimodipine. Regarding the ICA and BA arteries, this beneficial effect was even more pronounced. Although there was a tendency of better effects of the BAP-N group, regarding the overall effect in all territories combined, this failed to reach statistical evidence. In cerebral peripheral arteries, no differences were observed, and there was no difference in clinical outcome, too.

Identification of proteomic differences between squamous cell carcinoma of the lung and bronchial epithelium.

Proteins that exhibit different expression levels in normal and malignant lung cells are good candidate biomarkers to improve early diagnosis and intervention. We used a quantitative approach and compared the proteome of microdissected cells from normal human bronchial epithelium and squamous cell carcinoma tumors of histopathological grades G2 and G3. DIGE analysis and subsequent MS-based protein identification revealed that 32 non-redundant proteins were differentially regulated between the respective tissue types. These proteins are mainly involved in energy pathways, cell growth or maintenance mechanisms, protein metabolism, and the regulation of DNA and RNA metabolism. The expression of some of these proteins was analyzed by immunohistochemistry using tissue microarrays containing tissue specimen of 55 patients, including normal bronchial epithelium, squamous cell carcinomas, adenocarcinomas, and large cell carcinomas. The results of the immunohistochemical studies correlated with the proteome study data and revealed that particularly HSP47 and a group of cytokeratins (i.e. cytokeratins 6a, 16, and 17) are significantly co-regulated in squamous cell carcinoma. Furthermore cytokeratin 17 showed significantly higher abundance in G2 grade compared with G3 grade squamous cell carcinomas in both the gel-based and the immunohistochemical analysis. Therefore this protein might be used as a marker for stratification between different tumor grades.